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1. Introduction 

In 1966, Marc Kac asked his famous question, "Can one hear the shape of a drum?" [1]. 
This question can be rephrased as "does the Laplacian on every planar domain with 
Dirichlet boundary conditions have a unique spectrum?" . Ever since the time when Kac 
posed this fascinating question, physicists and mathematicians ahke have attacked the 
problem from various angles. Attempts were made both to reconstruct the shape of an 
object from its spectrum, and to find different objects that are isospectral, i.e., have the 
same spectrum. The interested reader can find an elaborate summary of these efforts 
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in PP-[I3|- In 1985, Sunada presented a theorem that describes a method to construct 
isospectral Riemannian manifolds [2] . Buser and later Berard expanded on this theorem, 
and offered a proof based on the concept of transplantation, as summarized by Brooks 
[21 m [5]. Over the years, several pairs of isospectral objects were found, but these were 
not planar domains, and therefore did not serve as an exact answer to Kac's question. 
In 1992, by applying an extension of Sunada's theorem, Gordon, Webb and Wolpert 
were able to finally answer Kac's question as it related to drums, presenting the first 
pair of isospectral two-dimensional planar domains P,[7j. Buser et al. later obtained 
a set of seventeen isospectral families of planar domains, both Neumann and Dirichlet 
isospectral [8]. Jakobson et al. and Levitin et al. extended the choice of boundary 
conditions by considering objects with alternating boundary conditions, and found sets 
of four planar domains that are mutually isospectral P, [10]. In the late 1990's, Gutkin 
and Smilansky reposed and answered Kac's question as it applies to quantum graphs 
[21]. Recently, Band et al. presented a pair of isospectral quantum graphs [26], whose 
construction was generalized to the method described in this paper. 

We begin by reviewing the terminology and the relevant definitions for quantum 
graphs. In section [3l we rederive the graphs constructed in [26], to help the reader gain 
an intuitive understanding of the method. Once the reader is familiar with the notions 
used, we formalize a theorem (section H]) along with a corollary, which together form 
the crux of the construction method. With the theorem in hand, we return to the basic 
example presented in section [3l and show that the isospectral pair can be expanded 
indefinitely — section [5l After describing the assembly process rigorously in section [6l 
we devote section [7] to further investigating the theoretical implications of the theory. 
Finally, in sections [HI [9] we demonstrate how to apply the construction to other types of 
objects, and present a variety of examples of graphs, drums, and manifolds. 



2. Quantum graphs 

A graph T consists of a finite set of vertices V = {vi} and a finite set E = {cj} of edges 
connecting the vertices. Each edge e can be identified with a pair of vertices {vi,Vk}- 
We denote by Ey the set of all edges incident to the vertex v. The degree (valency) of 
the vertex is = \Ey\. This becomes a metric graph if each edge is assigned a finite 
length /e > 0. It is then possible to identify an edge e with a finite segment [0, le] of 
the real line having the natural coordinate Xg along it. In this context, a function on 
the graph is a vector / = . . . , /le|£;|) functions f\^^ : [O, /gj ^ C on the edges. 

Notice that in general it is not required that for t> G ^ and e, e' G E^ the functions f\^ 
and /Ig, agree on v. 

To obtain a quantum graph, we consider the following Hilbert space: (F) = 
0i='i^^ ([0,/eJ) with the inner product: {f,g) = X]l='i /o"' f\e/9\ej dxe^. The operator 



which draws our interest is the negative Laplacian: — A / 

The domain of definition for this operator is the Sobolev space, W^'"^ (F), the space of 
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all functions / such that /|g G W'^''^ ([O, ZgJ) for all 1 < j < \E\. In addition we require 
the functions to obey certain boundary conditions stated a priori: for a vertex v E V, 
we consider homogeneous boundary conditions which involve the function's values and 



derivatives at the vertex, of the form ■ f 



complex matrices, / 



is the vector 







..,/ 






ei 







0. Here A^ and are x dv 

T 

of the vertex values of the 



(f ) is the vector 



function along each edge incident to v, and /' 

of outgoing derivatives of / taken at the vertex. Before stating the boundary conditions, 
the graph is merely a collection of independent edges with functions defined separately on 
each edge. The connectivity of the graph is manifested through the boundary conditions, 
which are local in nature: we relate the values of the function and its derivatives on each 
vertex, but no relation is assumed between those values on different vertices, or along 
the edges. In summary, a quantum graph is a metric graph equipped with a differential 
operator and with homogeneous differential boundary conditions at the vertices. One 
can generalize the metric Laplacian by including a potential or a magnetic flux defined on 
the edges. However, these generalizations will not be addressed here, and the interested 
reader is referred to the reviews [TU [15] . 

A standard choice of boundary condition which we adopt is the so called Neumann 
boundary conditioiij|: 



/ agrees on the vertices: Wv eV Ve, e' E : f 



f 



0. 



• The sum of outgoing derivatives at each vertex is zero: Vt; G V : ^ /' {v 

It is worth noting that a Neumann vertex of valency two can be added at (or 
removed from) any point along an edge without changing the eigenspaces of the 
Laplacian, and thus, from a spectral point of view, without really changing the graph. 
Thus, loops (edges connecting a vertex to itself) and parallel edges (edges with the same 
endpoints) can be eliminated by the introduction of such "dummy" vertices — we shall 
occasionally exploit this to simplify notation by assuming, without loss of generality, 
that we are dealing with graphs with no loops or parallel edges. A possible choice for 
matrices that correspond to the Neumann boundary conditions is: 

/l-l \ /00---0\ 



A,: 



\ 



1 -1 




B„ 



\ 






1 





1 / 



For a vertex of degree one, the Neumann boundary condition is expressed by the 
matrices A^ = (0) , B^ = (1) and means that the derivative of the function equals 
zero at that vertex. Another useful boundary condition for a degree one vertex is the 
Dirichlet boundary condition, which means that the function vanishes at that vertex: 
Ay = (1) , By = {0). It can be seen that the Neumann condition renders the Laplacian 



I This boundary condition is sometimes referred to as Kirchhojf condition in the hterature. 
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self-adjoint, which guarantees that its spectrum is reaL In general, Kostrykin and 
Schrader provide necessary and sufficient conditions that ensure the self-adjointness of 
the Laplacian: for every v ^ V the d^j x 2dv matrix {Ay By) must be of maximal rank dy, 
and the matrix Ay ■ must be self-adjoint [16]. For a quantum graph F we shall denote 
by $r (A) the space of complex functions on F (i.e., satisfying the boundary conditions 
at the vertices) which are eigenfunctions of F's (negative) Laplacian with eigenvalue A: 



which assigns to each eigenvalue its multiplicity. Two quantum graphs F and F' are said 
to be isospectral if their spectra coincide, that is ur = err'- 

Quantum graphs play an important role in the study of quantum chaos. This 
connection was ffist revealed by the work of Kottos and Smilansky [171 [IB]- They 
show that the spectral statistics of quantum graphs follow the predictions of random- 
matrix theory very closely. They propose a derivation of a trace formula for quantum 
graphs and point out its similarity to the famous Gutzwiller trace formula [IHl [20] for 
chaotic Hamiltonian systems. The trace formula for a quantum graph connects the 
spectrum of the graph's Laplacian to the total length of the graph and the lengths of its 
periodic orbits. The main result in the field of isospectrality of quantum graphs is that 
of Gutkin and Smilansky, [24j, where they use the trace formula to show that under 
certain conditions a quantum graph can be heard, meaning that it can be recovered 
from the spectrum of its Laplacian. The necessary conditions include the graph being 
simple and its edges having rationally independent lengths. When these conditions are 
not satisfied, isospectral quantum graphs indeed arise. An early example appears in 
[2T] , in which Roth obtains isospectrality exploiting a spectral trace formula. vonBelow 
[22] uses the connection between spectra of discrete graphs and spectra of equilateral 
quantum graphs to turn isospectral discrete graphs into isospectral quantum graphs. 
In [23] isospectrality of weighted discrete graphs provides isospectral quantum graphs 
whose edges vary in length. A wealth of examples is constructed in [211 125] , using an 
analogy of the isospectral drums obtained by Buser et al. [8j. A recent example is the 
pair of isospectral dihedral graphs presented in [26] . Their construction was generalized 
to obtain the more complete theory which is presented in this paper. 

3. A basic example 

Let F be the graph given in figure [Tj^a). The lengths of the edges are determined by the 
parameters a, 6, c and we impose Neumann boundary conditions at all vertices. G = D4, 
the dihedral group of the square, is the symmetry group of F. G consists of the identity, 
three rotations and four reflections. Let r denote the reflection of F along the horizontal 
axis and a the rotation of F counterclockwise by Y2. The axes of the reflection elements 



<l>riX) = {f eW'^T) I -A/ = A/} . 



(2.1) 



We define the spectrum of F to be the function 
(Tr : A 1-^ dime $r (A) , 



(2.2) 
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(a) 



(b) 



Figure 1: (a) A graph that obeys the dihedral symmetry of the square. The lengths of 
some edges are marked, (b) The same graph, showing the axes of the reflection elements 
in D4. 



in G are shown in figure [Dj^b). We can describe G and two of its subgroups Hi,H2 < G 
by: 

G = {o',t) = {e, a, a^, a^, r, rcr, rcr^, ra^} 
Hi = {T,Ta^) = {e, T, rcr^, a^} 
H2 = {ra.Ta^) = {e, ra, ra^, a^} 

Consider the following one dimensional representations of Hi,H2, respectively: 

Ri:{e^{l), r^(-l), ra' ^ {1) , a' ^ (-1) } (3.1) 

R2:{e^{l), ra^il), ra^ ^ {-1) , ^ {-1) } (3.2) 

We will use these representations to construct two graphs denoted by ^/ri, ^/r2 
(figure [2]) which will be found to be isospectral. 



N 



b 2a 



N -'c 

c 



2c 



D D 

(a) 

Figure 2: The isospectral graphs (a) ^/Ri, (b) ^/R2- Neumann boundary conditions are 
assumed if nothing else is specified. D stands for Dirichlet boundary conditions and N 
for Neumann. 

We now explain the process of building the quotient graph ^/ri. Let A G C and 
/ G $r(-^) be a function which transforms according to the representation Ri, i.e.: 



WgeHi, g~f = pR^{g)~f. 



(3.3) 
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In the l.h.s., the action of the group Hi on / is by Wx G F, {gf){x) = f{g~^x). 
We use the transformation law of / in order to deduce its properties: we know that 
rf = — /, which means that / is an anti-symmetric function with respect to the 
horizontal reflection. We deduce that / vanishes on the fixed points of r (marked with 
diamonds in figure [3]( a)). In a similar manner, we see that / is symmetric with respect 
to the vertical reflection, since ra^f = /, and therefore the derivative of / vanishes at 
the corresponding fixed points (the squares in figure [3]( a)). Furthermore, it is enough 
to know the restriction of / to the first quadrant (the bold subgraph in figure [3](a)) in 
order to deduce / on the whole graph, using the known action of the reflections: 

r/ = -/, Ta'f = f. (3.4) 

Our construction process is now complete. The quotient graph ^/ri is the subgraph 
which lies in the first quadrant, with the boundary conditions of Dirichlet and Neumann 
in the appropriate locations, as was found for / (figure [3|^b)). 




(a) (b) 



Figure 3: (a) The information we have on / which transforms according to Ri. Diamonds 
mark the vertices at which the function vanishes and squares the vertices with zero 
derivative, (b) The quotient graph ^/Ri which encodes this information. D stands for 
Dirichlet boundary conditions and N for Neumann. 



From this example, we conclude that the construction of the quotient graph is 
motivated by an encoding scheme. We choose a fundamental domain for the action 
of G on F, i.e., a minimal subgraph from which the entire graph can be reached by 
the action of the group. We take this domain to be the quotient graph ^/ri. We 
encode a function / G $r(A), which transforms according to the representation Ri, by 
a function / G $r/jj^(A). The encoding is described by the map ^ : $p^(A) $r/i,^(A), 
where ^^^(A) is the space of all functions / G $r(A) that transform according to the 
representation Ri, and is just the restriction map to the fundamental domain. An 
important observation is that given / G ^r/R^{X), it is possible to construct a unique 
function / G ^^^(A) (using fl3.4p ). whose restriction to the fundamental domain is / 
(this is the decoding process). It follows that is invertible and thus is an isomorphism: 

<l>^nA) = $rM(A) (3.5) 
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The quotient ^/r2 is constructed similarly. For any A G C, we consider / G $p^(A), 
which means that / transforms in the following way: 

Taf = f, ra^f = -f (3.6) 

and carry on with the same arguments as above. We do not specify the details of this 
process but rather summarize it in figure HI 




Figure 4: (a) The information we have on / G <I>p^(A). Diamonds mark the vertices at 
which the function vanishes and squares the vertices with zero derivative, (b) The quotient 
graph r/ij2 which encodes this information. D stands for Dirichlet boundary conditions 
and N for Neumann. 

The isospectrality of ^/ri and ^/r2 is a direct consequence of corollary 14.41 which 
appears in the next section. 

4. Representation theory and isospectrality 

Having informally exposed some of the relations between representations and quantum 
graphs, we begin a more precise examination, culminating in a general theorem on 
isospectrality Let F be a graph that obeys a certain finite symmetry group G; this 
means that the action of G preserves the lengths of the edges and the connectivity 
and boundary conditions at the vertices. We do not assume that G is the maximal 
symmetry group of F. For every A G C, $r (A), the vector space of all A-eigenfunctions 
of the Laplacian on F, is a carrier space of some representation of G. This follows 
from the Laplacian commuting with the symmetry group: V(7 G G, — A {gf) = 
g {— A f) = X{gf), which renders <l>r (A) closed under the action of G. Let R denote the 
representation carried by $r and decompose it into the irreducible representations 
of G. Such a decomposition allows us to present $r (A) as some direct sum of carrier 
spaces of the irreducible representations of G. Denote by Si, . . . , Sr the irreducible 

§ See [Appendix A| for a short review of the Algebra used in this section. 
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representations of G and assume that Si appears times in R, i.e., R = ^l^^UiSi. 
Then 

$r (A) = {vf' © . . . © Kt) ® • • • ® {^1^ © ... © v;^"-) 

r 

= 0<^'r'(A), (4.1) 

i=l 

where, for each i G {l..^'}, {Vf"^^^_^ are carrier spaces of the irreducible representation 
Si (one carrier space for each copy of Si in R), and their direct sum in $r (A) is denoted 
by (A). (A) is called the Si-isotypic component of $r (A); it is the vector space 
consisting of all A-eigenfunctions which transform under the action of G according to the 
representation Si. Denoting by 1g the trivial representation of G, (A) is called the 
trivial component of $r (A), and is the space of all A-eigenfunctions which are invariant 
under the action of G. 

We pause the algebraic discussion for the purpose of reexamining the example from 
the previous section. Recall that Ri is a representation of the group Hi, and that 
we have constructed a quotient graph ^/ri, such that the following isomorphisms were 
established: 

VA e C, ^ : <^^' (A) ^ <^r/n, (A) (4.2) 

By definition, the dimension of $r/jj^ (A) is crr/j,^(A), the multiplicity of A in the 
spectrum of ^/ri. In analogy, we denote by cxp^ (A) the dimension of (A). But what 
does this dimension tell us? It can be thought of as the spectrum that will be observed 
by someone who can only see functions which transform according to Ri. We call a^^ 
the Ri-spectrum of F, and note in particular that it is a subspectrum of T: a^^ < ar- 
In terms of dimensions (14. 2 p gives 

c^r' = • (4-3) 

We thus identify the role of the quotient graph ^/ri as having the same spectrum as the 
i?i-spectrum of the original graph. This will be the characterizing property of all our 
quotient graphs, and it is therefore time to generalize the discussion. 

As we have defined the i?i-spectrum by o"p^(A) = dim$p^(A), we now define the 
S-spectrum of F, for any irreducible representation S (not necessarily one dimensional), 

as 

which can be interpreted as the number of copies of S of which $p (A) consists. This 
is also equal to the number of copies of 5* in $r(A) (the corresponding in (14. ip ). 
Using the orthogonality relations of irreducible characters, we may rewrite (14. 4p as 
CTp (A) = (x5) X^rCA))^. We use this equality to generalize the definition: 

Definition 4.1. For any representation R, the i?-spectrum of F is 

f^r(A) := (xK,X$r(A)>G • (4.5) 
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Note that fl4.4l) holds for irreducible representations and is not true for all R (in 
fact, is even defined when R is reducible). Equipped with the notion of the 

i?-spectrum of F, we can now define what is a ^//? graph. 

Definition 4.2. A ^/n-graph is any quantum graph T' whose spectrum is equal to the 
i?-spectrum of F: 

^ — 

(Tr' = CTp . 

Since all r/_R-graphs have the same spectrum by definition, we allow ourselves, by 
abuse of language, to refer to the spectrum o/^/r. The following isospectrality theorem 
then follows: 

Theorem 4.3. Let T be a quantum graph equipped with an action of a group G, H a 
subgroup of G, and R a representation of H. Then ^/r is isospectral to r/indg/j. 

Proof. 

VA e C (Tr/niX) = cr^(A) = {xr , X^rW)H 

= (Xlndgi? ) X*r(A))G 

= [X) = CTr/indgijlAj 

where moving to the second line we have used the Frobenius Reciprocity Theorem. □ 

Corollary 4.4. If G acts on F and Hi,H2 are subgroups of G with corresponding 
representations Ri, R2 such that Ind^^-Ri = Ind^^-^a? then ^/ri and ^/r2 are isospectral. 

Remark. This corollary is in fact equivalent to the theorem, as can be seen by taking 
H2 = G,R2 = Ind|^i?i. 

The isospectrality of the pair of graphs constructed in section [3] (figure [2]) is obtained 
from the above corollary. Returning to that example, one first notices that the two 
graphs presented there obey definition 14. 2[ This is true since we have shown that 

VAgC $r/«^(A) ^ <l>f^(A) 

from which follows o"r/_R^ = o"p ^ and therefore the first graph presented in section [3] 
can be honestly called a V^i-graph. The same goes for the second graph, ^/r2. The 
representations R2 used in the construction satisfy the condition Ind^^i?i = Ind^2-^2 
and this enables us to apply corollary 14. 41 and conclude that ^/ri and ^/r2 are isospectral. 

5. Extending the basic example 

It is clear that theorem 14. 31 and corollary 14. 41 would yield isospectral examples only when 
the required quotient graphs indeed exist. We saw the existence of two such quotients in 
section [31 A proof of the existence of the quotient of any graph by any representation, 
along with a rigorous construction technique, are given in [27J. In this section, we 
present key examples which enable us to gain insight into this method, as well as to 
understand the procedure implemented in |27j . 
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We return to the basic example brought forth in section [31 wishing to extend it by 
discovering more quantum graphs which are isospectral to the pair of graphs ^/ri,^/r2 
in figure O Corollary 14. 41 offers a method by which this can be achieved: find a subgroup 
H3 < G and a representation i?3 of it such that Ind^^-Ri = lnd'^^R2 = Ind'^^R^. Then, 
^/rs is isospectral to ^/ri and ^/r2. Such a subgroup and a representation indeed exist: 

H3 = {e, a, a^ a^} 



'11 



cr 



-r 



a [—1] 



(5.1) 



3 _ V, 





B,u — 



(b) 



1 -i 





1 i 



Figure 5: (a) The information we haveon/ E (A). There is a factor of i between the 
values and the clockwise derivatives of / at the marked vertices vi,V4,. (b) The quotient 
graph ^/R:i which encodes this information. The only non-Neumann boundary condition 
is at vi and it is specified by the matrices Ay^ , By^ . 

We use the intuitive approach obtained from the basic example in order to construct 
the quotient ^/rs. Let / be a function that transforms according to the representation 
i?3. The action of the rotation element cr on / is given by 

af = tf. (5.2) 

This means that knowing the values of / on a quarter of the graph (for example, the 
quarter marked in bold in figure Et^a)) enables us to deduce the values of / on the 
whole of r. We therefore take this subgraph to be our quotient graph and check what 
boundary conditions we should impose on it. From (15. 2p we obtain 



/ 



^1 



^^4 



(5.3) 
(5.4) 



These equations suggest that we should identify vertices ■51,^4. They merge to become 
a single vertex vi in the quotient (figure EJ^b)). Equation (15. 3p gives the boundary 
condition for the values at this new vertex. The boundary condition for the derivatives 
at Vi is obtained from (15.41) . by recalling that / obeys Neumann boundary conditions 
at the vertex V4 of F: 



(5.5) 
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Therefore, the boundary conditions at Vi in ^/r-s are described by the matrices given in 
figure [5](b), and this concludes the construction of ^/r^. The isomorphisms 



VA e C, 



^ : 



(A) 



can be easily deduced from the construction process. 



(5.6) 

Taking dimensions gives 

Cr/j?3, which proves the validity of this quotient. 
We already obtained an isospectral triple, but this does not cause us to stagnate 
and we go further with our isospectral quest. By theorem 14. 3[ any ^/R-gTaph, where 
R = Ind^^i?!, would be isospectral to our three graphs. A simple calculation (See 



Appendix A) shows that R is the single two dimensional irreducible representation of 



G = D4. By a choice of basis we can describe i? as a matrix representation. Such a 
representation is: 



1 
1 



-1 
1 



1 

-1 



a 



ra 



1 

-1 



a 



ra 



-1 

1 



1 

1 



(5.7) 



We now construct the quotient ^/r. The graph obtained, which is shown in figure 
El^b) is the same as ^/ri (figure [2](a)) and we elaborate on this phenomenon later on. 
The construction process of this "new" graph is now presented for didactic reasons. 
The procedure is similar to those already described, although a slight complication 
arises from R not being one dimensional. 





N 
N 



•^1 Ji 



• D 

D 



(b) 



Figure 6: (a) Two copies of the graph T which assist in the visualization of the functions 
f^\f^^ on r. The fundamental domain is marked in bold, (b) The quotient graph ^/r 
formed as two copies of the fundamental domain. 



We consider two functions /{ , /2 € $r(A) (for some A G C) that transform 
according to the matrix representation (15.71) . It follows that f'^\ /g^^ G $p (A), and that 

and f^^ form a basis for a carrier space of the representation i?, which we denote 
by V^i^- We may proceed in this manner, choosing f^\f!f^^... such that each pair, 
I ff^ ) /i*'' I ) transforms according to (15.71) , and is linearly independent of the previous 



The Isospectral Fruits of Representation Theory: Quantum Graphs and Drums 



12 



ones. Therefore, each pair, |/i*\/2*j, spans a different carrier space of R, which we 

n 

denote by V/^ and we get that W — © K:'^- The number of carrier spaces is 



i=l 



n 



dim <l>^ (A) 
dim R 



(5.8) 
(5.9) 



(recall that R is irreducible). We wish to construct a quotient ^/r such that 
dim $r/ij(A) = n , 

which means that crr/fl(A) = ap (A). If this holds for every A, then ^/r is indeed the 
desired quotient graph (definition 14.21) . If we again relate the construction process to 
the encoding technique, we see that we can achieve (15.91) if each carrier space V/^ is 
encoded by a single function /^*^ G $r/ij(A), in a manner that {/^*^}"_-,^ is a basis for 

$r/fl. (A). We demonstrate this idea by encoding = Span < f[^\ f!^^ \. This is done by 



thinking of the basis functions f\^\ f!^^ as residing on two copies of the graph F (figure 
int^a)). Knowing the values of f^^ and f^^ on a fundamental domain for the action of 
G on r (e.g., the bold subgraphs in figure [n](a)) allows one to deduce the values of f^^ 
and f^^ on the whole graph, using the known action of the group (15.71) . Therefore, the 
quotient graph is the union of these two copies of the fundamental domain. Its boundary 
conditions can be concluded from (15. 7p . which gives the relations between the values of 



g^^^ between their derivatives: 



n 



(1) 



n 



i2 



(1) 



^1 



We recall that f\^^ satisfies Neumann boundary conditions at Vi. 



.^2 











n 



(1) 



[V4 



n 



(1) 



n 



(1) 



/ 



[V4 



[V4 



- n 



"(1) 



[V4 



Plugging (I5.12P into (I5.13p . (15.140 . gives the following relations: 



(1) 



^4 



n 



(1) 



^4 



(1) 



^4 



n 



(1) 



^4 



(5.10) 
(5.11) 
(5.12) 

(5.13) 
(5.14) 

(5.15) 
(5.16) 



Equations (I5.15P and (I5.16P motivate us to glue the two sub-graphs at the vertex and 
supplement the new graph with Neumann boundary conditions at this vertex, which 
we denote by ^4. The relations (15.101) . (15.111) give Neumann and Dirichlet boundary 
conditions, respectively, on the other vertices, and this fully describes the quotient 
graph r/ij (figure [6](b)). 

We now explain how this encoding enables us to prove that ar/u = dp , and ensure 
the validity of our ^/r. Given fi^\f2^^ as above, we form a function /*^^) e $r/_R(A) 



The Isospectral Fruits of Representation Theory: Quantum Graphs and Drums 13 

whose left half is equal to the restriction of f[^^ to the fundamental domain, and whose 
right half is the restriction of /g^"* to the fundamental domain. The considerations above 
apply for every carrier space V/^, and we can encode its basis /2*''| by a function 

(z $py^(A) in the same way as was done for i = 1. 

The set |J < fi \ \ forms a basis for $p (A), and is therefore linearly independent. 
i=i J 

We now show that {/*•*'' }^_-^ are linearly independent as well, which gives | dim $r (A) < 
dim$r/fl(A), i.e., dp (A) < ar/u^X). Assume that J2^=i'^if^^'^ = O5 so we have for 
the restrictions of | /i*"* , /2*'' | on the fundamental domain that ^"=1 Cj/[*^ = 

Sr=i ~ 0- Using the known action of the group (15.71) . which linearly relate the 

values of /2*''| everywhere on T to their values on the fundamental domain, we 

obtain Yl^^^i cJi^ = YJi=i = 0, hence, q = 0. 

In order to show the opposite inequality \ dim $r (A) > dim ^t/r (A) we employ the 
decoding process, which turns a function g on '^/r into a pair of functions {^1,^2} on V 
which transform according to (15. 7p . This decoding is done in the following way: first, 
we set cji to be equal to the restriction of g to the left half of ^Ir and ^2 to the 

FD _ _ FD 

restriction of g to its right half. Next, since (71,(72 should transform according to the 
matrix representation (15.71) . we know how to express the values of each of gi and ^2 on F 
as a linear combination of their values on the fundamental domain. We now show that 
starting from any set of linearly independent functions {(7*^*''}™ on ^/r, and performing 

the decoding process on each of them to obtain the set l^i*"*, ^2*"*! . of functions on F, 

the latter set is linearly independent as well. Assume that YlT=i (^i^9i^ + '^2^92^ = 0- 
Since | transform according to (15. 7p . we can apply the elements of G to this relation 

in order to find others; for example, from rcr^ we obtain Yl^i '^2^9i^ + c'i92^ ~ 0- 
Since R is irreducible, the matrices in (15. 7p additively span M2 (C) ([38], section 3C). 
Therefore, by applying a suitable combination of G's elements we obtain 

m m m m 

i=l i=l i=l j=l 

SO that YlT=i 9^"^^ = YlT=i ^2 ''fi'*'*'' = O5 therefore Vi c{^ = c^^ = 0. Hence 
idim$p(A) = dim$r/fl(A), which gives (Tr(A) = av/n^X), and finishes the proof (see 

(EH), dEl). 



As was already mentioned, the "new" graph ^/r we have obtained is in fact one of 
the graphs that we saw previously, namely, ^/ri (figure [2](a)). At first sight one might 
wonder whether theorem 14.31 has any practical use, as the isospectral quantum graphs we 
obtained from it are isometric. This fear is ungrounded. Notice that ^/r was constructed 
according to a specific matrix representation of i?, and we still have not checked how a 
change of basis (which changes the matrix representation) affects '^/r. There is indeed 
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a dependency on the choice of basis, as we now show. Another matrix representation of 
R is 




(5.17) 



We present only the matrices of the elements rcr^,rcr^ since we saw that they suffice 
to construct the quotient. Figure [7] is analogous to figure El for the current case. The 
fundamental domain that we choose is the same, and the difference will appear in the 
boundary conditions. 





(b) 




(c) 



Figure 7: (a) Two copies of the graph T which assist in the visualization of the functions 
/i, /2 on r. The fundamental domain is marked in bold, (b) The first stage in the gluing 
process. The encoding is done by restricting /i, /2 to the fundamental domain, (c) The 
quotient graph ^/r which is built using the matrix representation ()5.17p . The boundary 
conditions are described in (|5.22|) . (|5.23p . 



We consider two functions /i, /2 € $r(A) (for some A G C) that transform according to 
the matrix representation (I5.17p . The first column of the matrix representing rcr^ tells 
us that 



(5.18) 
(5.19) 



Evaluating fl5.18p . (15.190 on and using the knowledge that fi obeys Neumann 
boundary conditions on (see equations (15.131) . (I5.14p ) gives 



(l-v^/2) /i (^4) + 1/2/2 ,(^4) =0. 

(-l-v^/2) f[ {v,) + y2f, iv,)=0. 



(5.20) 
(5.21) 



This indicates how we should start gluing the two subgraphs in order to obtain the 
quotient. The first stage in this process, depicted in figure [7|(b), is to identify the vertex 
V4 in the two copies and turn it into the vertex of the quotient, with the boundary 
conditions that were derived in fl5.20p . (15.210 : 



A 




B 







1)4 



(5.22) 
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After treating vertices Vi, V2 similarly we obtain the quotient ^/r (figure [Tl^c)), whose 
remaining boundary conditions are given by: 

3/2 v^/2 




A 



B.„ 



B,, 




(5.23) 



This last example demonstrates that for a multidimensional representation the 
quotient graph depends on the explicit matrix representation by which it is constructed. 
One choice of basis for the matrix representation gave us a quotient which is identical 
to one obtained previously (figure EJ^a)), while another basis yielded a new quotient 
graph (figure [3(c)). As a matter of fact, all the examples of quotient graphs with 
respect to the various representations discussed so far (figures [2](a) and (b) and[5l^b)) can 
also be obtained as quotients with respect to the representation i?, by suitable choices 
of bases for its matrix representation. Furthermore, there are many other quantum 
graphs isospectral to these. For example, we consider an arbitrary orthogonal matrix 
representation of i?, which is parameterized by : 



TO" 



cos^ 9 - sin^ 9 
— 2 cos sin 6* 



2 cos 9 sin 9 
cos^ 9 - sin^ 9 



—2 cos 9 sin 9 
- cos^ 9 + sin^ 9 

cos^ 9 — sin^ 9 
— 2 cos 6' sin 



(e.g., (15.171) is obtained from 9 = ^/s). Using the general construction method, which 
is described in the next section, we obtain from this matrix representation the quotient 
given in figure [7](c), with the following boundary conditions: 



A,. 



( 2sin2^ 
sin 20 








5,, 








V 



2-2 sin^ 9 
-sin 20 



sin 29 \ 
2-2 sin^ 9 






- sin 29 

2 



A,. 



( 1 - sin 29 
2 sin^ 0-1 




2 sin^ - 1 \ 
1 + sin 20 




1)4 



2sin^0 j 

















v 



1 + sin 20 
1 - 2 sin^ 



1-2 sin^ 
1 - sin 20 j 



The matrices above are not square ones as is required from the definition of the boundary 
conditions. However, since their role is to describe linear restrictions and due to the fact 
that they are all of rank one, they can be reduced to square matrices by deleting the 
appropriate row^. We thus get a continuous family of isospectral graphs. We already 
met some of its members. For example, = gives the following boundary conditions: 



A 



A 



V2 



A 



V4, 




(5.24) 



See also the discussion that comes after (16.11). (16.21). 
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B.. 





1 1 



) 



5, 





2 



) 



(5.25) 



When applying this to figure [Tl^c), we notice that does not remain a vertex of degree 
two, but rather, sphts into two vertices of degree one, one having Dirichlet boundary 
condition and the other Neumann. The vertices vi and f2, however, stay connected and 
obtain Neumann boundary conditions. The resulting quotient is thus the one that 
we have already obtained as ^/r2 (figure [2](b)). In a similar manner, the quotient 
(figure [2](a)) is obtained from the choice 6 = We conclude by pointing out 
that the graph described in figure [7](c) is a good prototype for the isospectral family 
mentioned, yet it also might be misleading, since some members of the family have 
boundary conditions that tear apart the edges connected to some of the vertices, and 
thus change the connectivity of the graph. One should also pay attention to the fact 
that we have treated only orthogonal representations of R. We may further extend 
the isospectral family presented above by considering the broader case of all matrix 
representations of R. In particular, the quotient ^/r-j, (figure [5](b)) is obtained from the 
unitary representation 



6. The rigorous construction of a quotient graph 

The rigorous formalism of the quotient graph construction is fully described and 
proven in [2?]. Here we summarize this method in accordance with the discussion 
and demonstrations presented so far in this paper. Let F be a quantum graph with a 
finite set of vertices V and a finite set E of edges connecting the vertices. Let G be 
a finite group that acts on F, and R a (i-dimensional representation of G. We assume 
for now that G acts freely on the edges, i.e., ge ^ e for e E E, id g E G, and 
leave the treatment of a non-free action on the edges for later (section 17.21) . We may 
choose an ordered basis B = (fej)^^]^ for R, with respect to which we think of it as a 
matrix representation. We choose representatives {e*}[^^ for the orbits ^/g, and likewise 
{vk}t=i for v/g. We shall assume, by adding "dummy" vertices (vertices of valency two 
with Neumann boundary conditions) if needed, that G does not carry any vertex in 
V to one of its neighbors. This ensures that no edge is transformed onto itself in the 
opposite direction, which would force us to take half the edge as a representative. The 
quotient graph ^/r obtained from these choices is defined to have {ffcl^Li &s its set of 
vertices, and {ej}*._^^"^ for edges, where each e* is of length l^i. The vertices Vk, Vk' of 
^/r are connected by the edge e* if there exist g,g' E G such that e* connects gvk to 
g'vk' in F. In such a case the vertices Vk, Vk' are connected by all the edges {e*}^_^. 
Until now we have used only the action of the group G on F in order to determine 
the edges and vertices of the quotient graph and its connectivity. Now we also need to 
use the information from R, B, and the boundary conditions of F, in order to specify 
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the boundary conditions at each vertex in ^/r. Let the boundary conditions at Vk 
in r be described by A^j^, B^^. From here to the end of this section we focus on the 
vertex Vk and its predecessor Vk, and keep in mind that our parameters depend on k 
even when it is not reflected by the notations. The set of edges incident to the vertex 
Vk can be written as Eyi_ = {gie'^'-}^=i, for some {gi}^^i in G. Note that repetitions can 
occur among the z/;'s, and also among the ^^^'s, and that their total number is n = da^. 

The set of edges entering Vk is E^i^ = je^'' } !<;/<„ , where {yUpjllLi is defined to be the 

i<j<d 

set of distinct values among {i^iYi^i- Obviously, m < n, and the relation between the 
sets {fJ,l'}]?^-^^, is given by the n x m matrix 

e' , = / ^ = /^c 

I otherwise 

The resulting matrices (see |27j), describing the boundary conditions at Vk, are: 

A,^ = {A,^ ® h) ■ diag {[pR {gi')]^, ...,[pr {g-')]Bf " (©' ® Id) (6.1) 

5,, = (5,, ® /,) ■ diag ([pH {gi')]s^ ■■■^[PR {9n')]Bf ■ (©' ® Id) ■ (6-2) 

Observe that the boundary conditions (16. ip . (16. 2p encapsulate the information of the 
boundary conditions on T (given by A^j^^, B^J, the action of the group (6' and {(?«}), the 
representation (pr) and the basis that was chosen for the representation [B). However, 
they may fail to be square matrices as our definition of a quantum graph calls for. We 
therefore need to discuss the question of how many linearly independent restrictions 
are dictated by those boundary conditions, i.e., find rank (^1^,^. S^^). Let us assume 
that the original boundary conditions on Vk were linearly independent, i.e. (Ay^Bui^) 
is of maximal rank, (i^j,. If the action of G is free not only on the edges but also on 
the vertices, then 0' is a permutation matrix and we get that At,^. and B^^ are square 
matrices. Furthermore, rank {A^^, B^^) = d ■ rank (A^j. By^) = d ■ dj,^ = dy^,, which means 
that {Ayi^ By^) is also of maximal rank. 

In the general case, it is shown in [27j that if the matrix Ay^, ■ B\^ is self-adjoint 
then ia,'nk{Ay^ By^) = dy^. This means that we may eliminate rows from [Ay^By^), 
and remain with square Ay^, By^ which still describe the same boundary conditions. 
Therefore, starting from a graph F whose Laplacian is self-adjoint, any ^/r produced by 
this construction method would be a valid quantum graph (whether the Laplacian on 
^/r is also self-adjoint is discussed in section 17731) . This also demonstrates the benefits of 
using quantum graphs to implement the quotient construction. Although one may apply 
the above procedure to other types of objects (see section [9]), the resulting quotients 
might not be objects of the same type as their predecessors. The advantage of quantum 
graphs in this context is that the self-adjointness condition guarantees that the quotient 
is a quantum graph as well. 

We now recall the intimate relation between the construction of ^/r and the 
encoding process, which for each A G C takes functions /i, . . . , G $r(A) that transform 
according to R, and maps them to a single function / G $r/ji(A). This encoding is 
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described by: 

(6-3) 

J 

This relation gives the decoding process as well — starting from / G ^r/ii{X), one 
reconstructs the values of fi, . . . , fa G $r(A) on the edges' representatives {e*}. The 
values on the rest of the edges are then determined by the action of G and the 
representation R. 

We elaborate on the functions fi,---,fd&^rW that play a role in the decoding- 
encoding process. For two representations R, R' with corresponding carrier spaces V, V , 
Home (V, V) is the vector space of all linear transformations T : V ^ V that respect 
the action of the group G, i.e., 'ig E G,v E V : T{pR{g)v) = pRi{g)T{v). Such 
linear transformations are called intertwiners of the representations i?, R! . It is known 
that {xr-iX^tw) G ~ dimHomc (i?, $r (A)) ([38], section 9C) and we recall that the 
quotient graph's construction guarantees av/R^X) = (xR; X^rCA))^- We therefore have 
that $r/fl(A) = Home (i?, $r (A)), which offers another point of view on our quotient 
graphs. We demonstrate this by an example. Let G be a group and R = S Q) S Q) S, for 
some irreducible representation S of G, and let A G C be such that $r(A) = Vf © V2 , 
where Vf, V2 are carrier spaces of S. Since dim Home (i?, $r (A)) = (xj?) X<i>r{A))(j = 
(3xs, '^Xs)g ~ 6' ^^"^ choose six linearly independent intertwiners from to •I'r (A). 
One such choice of intertwiners can be described as follows. We decompose a carrier 
space of i?, V^, into three carrier spaces of S and each intertwiner then sends one of 
these three copies onto either V-^ or V2 ■, and the other two copies to zero. 
The quotient '^/r is constructed with respect to a certain basis B of R. The 
image of this basis by each of the above six intertwiners is a set oi d = dimi? 
functions, /i, . . . , G $r(A) which transform according to R 0, and this set can be 
encoded by a function / G $r/j,(A) according to (16.31) . We thus obtain six linearly 
independent eigenf unctions of the eigenvalue A on the quotient '^/r, which demonstrates 
the equivalence $r/fl(A) = Home (-R, $r (A)). 

7. Application of the method — Further investigation 

Having revealed the key elements of the construction method, we are ready to present 
its theoretical implications, as well as various issues that may interest those who wish 
to apply it. 

7.1. Sunada's isospectral theorem 

In his famous paper [2], Sunada provides a method for producing isospectral Riemannian 
manifolds. Phrasing his result somewhat loosely, given a manifold M equipped with a 
free action of a finite group G, and subgroups Hi, H2 of G which are almost conjugate, 
the manifolds *Yhi and *Y//2 are isospectral. Hi and H2 are said to be almost conjugate 

% Note that in this example, this set of d functions is not a hnearly independent set. 
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if V(7 G G \{[g]r\Hi}\ = \{[g]r\H2}\, where [g] denotes the conjugacy class of in G. As 
mentioned and proved in [5], i^i and H2 are almost conjugate iff Ind^^l//^ = Ind'^^lnzi 
where Ijj. is the trivial representation of Hi. Comparing this to corollary 14.41 we see that 
Sunada's theorem is obtained as a particular case of the corollary if ^'^/Hi is isospectral 
to ^/'i-Hi- This is indeed the case since the space of functions on ^Y^i is isomorphic to 
the space of functions on M which are invariant under the action of Hi. But the latter is 
the trivial component of the space of functions on M, which by definition is isomorphic 
to the space of functions on ^Viff») so that 

as claimed. During the discussion above we have applied our isospectral theory to 
manifolds, as opposed to quantum graphs, for which it was developed in the previous 
sections. This should not bother us - the method can be applied to other geometric 
objects, as will be demonstrated in section [9l Finally, we note that the equivalence of 
Ind^^l//^ = Ind^2li?2 Hi, H2 being almost conjugate was already used by Pesce [28] 
to give another proof for Sunada's theorem. Our application of Frobenius Reciprocity 
for the proof of theorem 14.31 is similar to that of Pesce, whose proof is also summarized 
comprehensibly in [5]. 

7.2. A non-free action on the edges 

We now treat the case of a group whose action on the edges is not free. We have already 
encountered examples of non-free action on the vertices, and saw that this may result in 
the vertex in the quotient graph having a degree smaller than that of its predecessor. An 
action which is not free on the edges results in some more interesting features. Denote 
by F the quantum graph which is an equilateral tetrahedron equipped with Neumann 
boundary conditions at all its vertices (figure [8]) . 




Figure 8: An equilateral tetrahedron with Neumann boundary conditions. It is equipped 
with an action of S4, permuting its vertices vi, V2, V3, V4. 
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Note that we have added "dummy" vertices in the middle of the edges, as described in 
the previous section. acts by permuting the vertices Vi, V2, V3, v^, and this determines 
its action on the whole of F. Let S denote the sign representation of S'4: 

WaeS, ps{a) = {sgna) . (7.1) 

We choose the edge and the vertices Vi, as a fundamental domain for the action of 



S'4. A function / which transforms according to S satisfies (3 4) / = ps [(3 4' 



Evaluating this on the edge e gives / 



/ 



(3 4)el 



(3 4)/ 



-/ 



/ = -/■ 
, which 



means that / vanishes on and therefore also on the whole of F. We conclude that F's 
Laplacian has no eigenfunctions which transform according to S, so that ^/s is just the 
empty graph. Putting this another way, the edge of F was chosen as a representative 
for the orbits of the edges under the action of the group, ^/s4. We therefore expected 
to have a single edge e\ in the quotient. However, this edge has disappeared due to 
the action of the stabilizer of e^. Specifically, the stabilizer is Ggi = {id, {3 4)} and 
the restriction of the representation on it is not trivial: ps 7^ 1- This caused the 

disappearance. 

Moving to a more complex demonstration of this principle, we examine the 
permutation representation of S'4, which we denote by R. Some of the matrices of 
R using its standard basis are: 



(7.2) 





(2 3 4) 




(2 4 3) 




We choose the same fundamental domain and consider the values of the functions 
/i 5/2 5/3 5/4 (which transform according to (17. 2p ) on it. We now have four edges, 
e\, el, 63, e\, that are supposed to form the quotient graph after we properly glue them 
and deduce the boundary conditions. From the action of (2 3 4) and (2 4 3) on /i, 
the chosen matrix representation for R, and the fact that vi has Neumann boundary 

[vi) = 0. Likewise, from the action of (1 2) on /i we 



conditions, it follows that f[ 



conclude that /i 



/2 



(fs) and f[ (fs) + (fs) = (where the derivatives 



are outgoing from v^). This allows us glue the edges e\, el and deduce a Neumann 
condition on the left side of e\ (figure [9]). These considerations are similar to those 
used before, and the new part comes when we observe that (3 4)/3 = f^. This gives 



/a 

64) in 



/4 



and means that we do not need both of the corresponding edges (i.e., 63 and 



the quotient graph El- We notice that the problem arises since the representation 
restricted to the stabilizer Ggi = {id, (3 4)} has a non trivial component (p^ [(3 4)] has 



This is because the quotient construction is motivated by encoding the functions on T which transform 
according to R, without redundancies. 
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an eigenvalue not equal to one). The solution for this situation is to change our basis 
in the following way: 

/i = fi, /s = 75 ( /s + /4 
/2 = f2, /4 = "^75 ( /a ~ fi I , 
This forces us to rewrite (17. 2p as 



:i2) 



(2 3 4) 



10 
10 
10/' 
1 



10 

i/V2 

yv2 1/2 

l/^/2 -1/2 





-1/v^ 

1/2 
-1/2 



(3 4) 



(2 4 3) 



10 
10 
10 

-1 



10 

1/v^ Vv^ 

1/v^ 1/2 -1/2 

-1/^2 1/2 -1/2 



(7.3) 



Examining the action of (3 4) on /4, we get 



and this turns out to be the edge 



and /2 



are 



which disappears. The previous considerations we made concerning fi 

still valid since /i = /i, /2 = f2- The last step in the gluing process follows from the 
second column of the matrix representations of (2 3 4), (2 4 3), and the evaluation of the 
corresponding relations on e^: 

h 
h 





+ VV2/4 

el 


el 


(7.4) 


) = 


- VV2/4 

fii 


el 


(7.5) 



Recalling that /4 



(234) /2 
(2 4 3)/2^ 

and using (17.41) . (17.51) to obtain relations for the derivatives as 



well, we finish the process and give the quotient r/5 that is shown in figured 
1 1 _i A 

B 



J- 111 



N - 



- N 



V2 -1 






1 



Figure 9: The quotient ^/s. The only non-Neumann boundary condition is at vi and it 
is specified by the matrices Ay^ , B^-^ . 



We now return to the general construction procedure described in the previous 
section, and discuss the adjustments that should be made in order to apply it to the 
case of a non-free action on the edges. After choosing the representatives {e*}[^j^ for the 
orbits ^/g, we consider for each i the representation Res^ .S* with its carrier space V^. 

We decompose into its trivial component, (V^) ^ and its complement, (V^) . 

e'' 

We denote the dimension of the trivial component of by di and choose a basis for 
it: {fejl^.^-^. We complete this to a basis for the whole of by adding vectors from 

the complement of (V*) "^e^ and we denote this basis, \b]} ._ , by B\ Let < /.• \ 

J ' J— 1 I eO j=l 
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be functions on the edge e* which transform according to the representation R as given 
by the basis i?*. As in the preceding examples we have that Vj > c/j, /■ =0 0- 
Therefore, in the assembly of the quotient, we endow it with only di copies of the edge 
— the copies which correspond to the functions < \ . Equations (16.11) . (16. 2p are 

l eO j=l 

replaced in the general case by the following: 

A, = {A,^ ® Id) ■ diag ([pfi {9i')]l., ,---,[pR {9n')]l.^y ■ (7.6) 

B^^ = (5,, ® /,) ■ diag ([p^ {9l')]l., ,...,[PR {9n')]l..y ■ (7.7) 

One obvious difference is the use of a separate basis for each edge that enters the 
vertex. The other change is due to G which is defined to be the nd x dy^^ matrix obtained 
by removing from (0' ® Id) the columns {{i ~ 1) ■ d + j} i<i<m ; these are the columns 

which represent the functions whose vanishing on the corresponding edge (e*) of the 
graph r causes the disappearance of (d — di) copies of that edge from the quotient. Note 
that (16.11) . (16. 2p are obtained as a special case of (17.61) . (17.71) when the action is free on 
the edges and a single choice of basis is made (i.e., B^ = . . . = B^ = B) . In addition, the 
encoding\decoding process in this case is then slightly altered from the one described 
by (16.31) . For each A G C and 1 < z < /, a set of functions /{,... G $r(A) that 
transform according to the basis B^ of R is mapped into a single function / G $r/fl.(A) 
by: 



yl<J<d^, f 



Note that the functions < /, 



l<i<7 
d^<j<d 



(7.8) 

are zero and therefore not encoded. Finally, the 



discussion that comes after equations (16.11) . (16.21) regarding the rank of A^^^B^^ is valid 
here as well: whenever the Laplacian on F is self-adjoint, the produced quotient obeys 
the definition of a quantum graph. 

7.3. The self-adjointness of the Laplacian on the quotient graph 

A natural question to ask is whether the Laplacian on the quotient graph produced by 
the above construction is self-adjoint. The necessary and sufficient conditions for the 
self-adjointness of a quantum graph were described in section [2l Their examination in 
light of (17. 6p . (17. 7p gives the following: 

Proposition 7.1. Let ^/r be a quotient quantum graph constructed as explained in 



section {7.^ . 

(i) If T's Laplacian is self-adjoint, G acts freely on F (both on the edges and the 
vertices), and [pr {g)]% is unitary for all g E G, I < i < I , then the Laplacian on 
r/_R is self-adjoint. 

* This can be deduced from observation (|7.10p . 
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(a) If T has Neumann boundary conditions and [pRig)]^ unitary for all g,i, then 
the Laplacian on ^/r is self-adjoint if and only if for every vertex v at least one of 
the following holds: 



(a) 



0. 



(b) All stabilizers {Gi}-^^_ are of equal order. 

Remark. Since G is assumed to be finite, one dimensional representations are unitary 
in all bases. 



Proof. We start by recalling the observations made in sections El I7.2[ When T has a 
self-adjoint Laplacian, our construction ensures that for every vertex v of the quotient 
we have rank (^4^, B^) = d^. Therefore, one condition required for the self-adjointness of 
the quotient is fulfilled, and we are left to verify the other, • Bl = ■ A\j. 

(i) Let V he a. vertex in F, with E;^ = {gie^'}^^i, and v the corresponding vertex in 
^/r. Since the action is free both on the edges and on the vertices we have 0' = In 
(by appropriate indexing of the edges) and 6 = 6' ® /d = Id^ (since the free action 
also implies = dj, ■ d = n ■ d). The unitarity of the representation with respect 



to the bases 5* and B gives that D : 
a unitary matrix. Therefore, 

A^-Bl= {A^ (g) Id 

= iAi,0ld 



diag (^[pr {g^ , • • • , [pij (^n^)]?-™) is 



B,~, ■ Al 



(fic ® 
^Id = 



a _ 



Id)-D- 1, 



dv J 



idV = {a,-bI 

-B,, ■ a], . 



Let V be a vertex of ^/r, and v its predecessor in F. Recall (section [6]) that 
Ev = {fi'iC'^'liLi' ^hat attain in total m distinct values, {pi}^i- We 

assume (by reordering, if necessary) that {z/j}™^^ are different: this means that 
{gi(i^'}T=i representatives for the different orbits of the edges in For each 
1 < i < m we choose a set of representatives {tP^y^_^ (rj = [G^, : Gg^e-'i]) for the 
left cosets of Gg^e^i in G^u, so that when writing E^j = [tlgiC^^^ i=i..m we obtain each 

j = l..rj 

edge in E^ exactly once. This is the ordering with which we represent the boundary 

conditions at v. 

Denoting 



0ij 
6 = 



, -1 



PR {^Idi) 
diag (01,1, . . . ,01,^^,02,1, ,0r 















G Mrfxd (C) 
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we have 



91,1^'isLi 



91,1^1 sl.n 



gm.lPm 0m. 1 



0m,lPm0m.?> 



The Neumann boundary conditions at v can be expressed by 



and using; (USD, (1Z21) we have B^Al = {B^ ® Q (2500t(St (A- o /^)t and can write 
S^Aj, in d X d blocks: 



1 


-1 






^° 


■ 


















1 








• 
















1 ■ 





B.Al 

[Mi,i~Mi,2] ■■■ [Mi,i~Mi,r,] [Afi.i ~ Af2,i] 

where for 1 < i < m, 1 < j < r j we denote 



fc=i 

Due to the fact that B^jAl is strictly lower-triangular, it can only be self-adjoint 
if it is zero, i.e., if Mj j is the same for all i,j, and we wish to analyze when this 
occurs. First, we note that 

by the following observations: 

(a) For every irreducible representation of a group G we have 



1 ^ / / S^l 

' ' 9<^G I ^ 



(7.10) 



(b) As explained in section 17.21 we chose B"^ = {b'^'}^_^ such that 

is a basis for the trivial component of Res^,^ .R, and is a basis for 

(Res^^^ 



By our unitarity assumption 



PR {tl9i) 



1 B 



[pR {^i9i)]B ' SO that 



k=l 



B 
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Using (17.91) we find 



k=l 



B 



k=l gdG^Vi 



and since giGg'^ig, 



-1 



|G 



-1 



fc=l geG, 



\Gz"i 



E ^9)] 



B 

B ' 



as for every z we liave tlie disjoint union G^ = IJfc=i tiGgiii^, , and tj G for all i 
and j. We therefore have Mij = Mi^j' iff either {Gg'^i \ = IGg"*' | or J2geGy PR (9) ~ ^■ 
Observation (17.101) shows that the latter happens if and only if Res^.i? has no trivial 
component, i.e., {xr, 1)^- = 0. 



□ 



Corollary 7.2. IfT 's Laplacian is self adjoint, and G acts freely on T, then a quotient 
with a self-adjoint Laplacian can be constructed for any representation of G. 

Proof. This follows from the first part of the proposition; as the action is free, we have 
no restrictions on the bases {-B*}^^^ used in the construction. By Weyl's unitary trick, 
any finite dimensional representation of a finite group is unitarisable (see for example 
[37] ) and therefore we can choose a single basis B with respect to which i? is a unitary 
matrix representation, and take = . . . = = B. □ 



7.4- Transplantation 

The transplantation technique is one of the simple and illuminating proofs for 
isospectrality. For two isospectral objects, the transplantation is a linear map which 
bijects each eigenspace of the first object's Laplacian onto the corresponding eigenspace 
of the second. The existence of a transplantation proves that the spectra of the two 
objects coincide, including multiplicity. One usually encounters the transplantation 
technique when deahng with isospectral objects consisting of building blocks. Several 
copies of a certain fundamental building block are pasted together along their boundaries 
to form the complete object. The second isospectral object is composed of the same 
building blocks, pasted in a different manner. The transplantation then presents the 
restriction of an eigenfunction to a building block of one object as a linear combination 
of such restrictions to several building blocks of the other object. The linearity of the 
eigenvalue problem ensures that the transplanted function is indeed an eigenfunction 
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of the second shape, as long as it obeys the boundary conditions. The idea of 
transplantation is further discussed and demonstrated in [31 HI O [8], [9], [TOl [291 130 ]. 

In order to discuss the notion of transplantation as it arises from our method, it 
is crucial to examine the exact procedure by which the isospectral quantum graphs 
are constructed. The main tools for producing isospectral graphs are theorem 14.31 and 
corollary 14.41 However, we have already noticed that ^/r is not a single graph, but rather 
a continuum of isospectral graphs, each constructed by a certain choice of basis for 
R. Therefore, isospectrality emerges also from a change of basis for the representation. 
Furthermore, even when one works with a single basis, ^/r might depend on the choice of 
representatives for the orbits ^/g, ^/g. An example of producing isospectral graphs using 
this freedom is shown in section [8l2l This variety of degrees of freedom (representations, 
bases, representatives) motivates us to address each of them separately and one can 
combine them appropriately in order to construct the over-all transplantation. 

Proposition 7.3. Let T be a quantum graph equipped with an action of a group G, and 
let R be a representation of G. Then any two ^/r graphs constructed according to the 
recipe given in sectionslBi \ 7.S\ are transplantable. 



Proof. We deal with the following degrees of freedom: 

(i) The freedom to choose a basis. 

Let {e*}^^^ be representatives for the orbits ^/g of T (the choices of representatives 
for the vertices do not affect the transplantation). For every 1 < i < J, let B^, 03' 
be two bases for R chosen to correspond to the edge e* (following the prescription 
given in the end of section 17.21) . Let Fi be the quotient obtained as ^/r when 
we choose the bases {B^}^^^ for R and F2 the quotient ^/r obtained by the bases 
{Q3*}f^]^. Our motivation is the following: for every A G C, $ri(A) and $r2('^) serve 
to encode the same space (if R is irreducible then this is the i?-isotypic component of 
$r (A). In the general case, it the corresponding space of intertwiners) . Therefore, 
a bijection between them can be obtained by composing the 'Ti-decoding" with the 
'T2-encoding" . However, since functions on the two quotients encode functions on 
F which transform according to different bases, a suitable change of basis is required 
between the decoding and encoding. Starting with / G $ri(A), and denoting the 
corresponding function in $r2('^) by ip, we obtain 



/ 



where [I]^ is the change of basis matrix and \ fj\ , {^k}'k=i functions in 



k=l 

t \mS. 



k=l 



-I k,j 



,^k 



-I k,j 
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$r (A) described by (17. 8p . By the construction, transform according to B\ 

and {lpk}k=i according to 03*. We now note that the bijection we have constructed 
is indeed a transplantation, as it represents the restrictions of one function to each 
edge as a hnear combination of restrictions of the second function to edges of the 
same length. We also note that the choice of the global basis B (see (17. 6p . (17.71) ) 
does not affect the transplantation. 

(ii) The freedom to choose representatives. 

Let {e*}^^^ , be two sets of representatives for the orbits ^/g. We describe 

the connection between them by choosing such that g^e^ = for 1 < i < I. 
For each i, let = be a basis for R chosen to correspond to the edge 

e\ A natural choice for a basis 03* = which corresponds to the edge 

would then be = g^b^ Let Fi be the quotient obtained as ^/r by the choice of 
the representatives {e'^}l^i and the bases {B^}1^^, and F2 the quotient ^/r obtained 
from and {^^jf^^. Working with these bases for F2 does not limit generahty, 

as other choices can be handled by the first part of the proposition. The merit of 
{Q3*}f^]^ is that it gives an extremely simple transplantation between / G $ri(A) 
and the corresponding ip G $r2('^)- 



/ 



g'f] 



■ - =^3 



We have used (I7.8P on the first and last equalities. The second follows from the 
action of G on $r (A), and the third from the choice of bases for the matrix 

representation (as before, . ; {"^jlj-i are functions in $r (A) that transform 

according to the matrix representations of R as given by the bases i?*, 03*). 

□ 

Proposition 7.4. Let T he a quantum graph equipped with an action of a group G, 
H a subgroup of G, and R a representation of H. Then there exists a transplantation 
between every two graphs ^/r and ^/indgfl which are constructed by the recipe given in 
sectionsl^ \ 7.2. 



Proof. The outline of the proof is similar to that of proposition 17.31 We start by 
describing a convenient choice of representatives and bases for the construction of each 
of the quotients ^/r and ^/ind%R. For every A G C we obtain a connection between two 
sets of functions in $r(A), which transform according to those two sets of bases. This 
connection yields the transplantation between the functions in $r/j;(A) and $r/indgii(A)- 
Let {e'^Yi^i be the representatives for the orbits ^/g oi T used in the construction of 
r/indgfl. We deal separately with each representative £* and address the simpler case 
of Ggt = {id} first. Let -B* = {&*}^._-^ be a basis for the representation R. We choose 

representatives for the left cosets of H in G: {tfcjjj'lf . A possible basis for Ind^i? 
tt See I Appendix A| 
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is 03* = \tkb]} i<k<[GH]- Let A G C. For < f^> , a set of functions in $r(A) which 

^ ^ i<i<d I J l<i<d 



transform according to B\ we consider |v^(fcj) = i<fc<[G./f] ' ^ set of functions in 



i<j<d 



$r(A) which transforms according to QSj. We handle the case in which ^/r is constructed 
by choosing the representatives for the orbits '^^'/h of T to be je*^*''^-' = t^^e'^y^_^^\ for 
which we obtain the transplantation: 



=^M.=^^5.=5.-.. = 5.,,,=/ (7-11) 



We now treat the case of a non-trivial G^i. Obviously, the construction of each 
quotient depends on the various choices for bases and the edge representatives. We 
proceed in the following order: starting with any representatives for ^/g, we pick 
representatives for ^/h, {(i(i^k)}- We then assume we are given any bases for R, {B^''''^}, 
which fit {e(i,fc)} and use them to produce bases {23*}, which fit {ii}. The instrument 
which enables us to advance in this manner is the double coset structure. 
We start by choosing representatives for the {G^i,H) double cosets in G, {tkYi^^i- 
then denote {t(fc,i)}'^_-^ = {^A:}fc=i and complete this to a set of representatives of the 
left cosets of H in G, {t(k,m)} i<k<i {Ylk=i^k = [G : H]), obeying the condition 

\/l<k<l, Vl<mi,m2<nfc, 3geGg^ s.t. t(^k,mi) = gt(k,m2)- 

Note that the choice of {t(^k,m)} obviously depends on e\ but we omit the i index here 
to simplify the notation. Returning to the graph, the condition above means that 

V1<A:</, Vl<mi,m2<n,, t-,]^^f^ = t^,]^^f\ (7.12) 

It is now possible to describe the representatives for the orbits ^/h, used in r/^'s 
construction. For each 1 < i < / we can choose the representatives for the orbits '^^'/h 
as = {^fc^^*}i<jfc<r denote these representatives by e'^*''^^ = t^-^e^ 

The union of all these, je*^*''^^} i<i</ , forms our choice of representatives for the orbits 

i<fc<i 

e/h. 

We now assume that for each k we are given i^^*'*^) = , a basis for R which 

corresponds to the edge e^^'^\ Namely, is a basis for the trivial component 

of Res^ R, and Ifel*''^^! is a basis for ffi ^Res^^ R] (note that du depends 

on i, but we do not mention this to simplify the notation). A possible basis for Ind^i? 
is \ t(^k,m.)b^j''''^ \ i<m*<4 5 but the construction of ^/ind%R requires choosing a basis which 

' l<j<d 

corresponds to the edge in the sense mentioned above. By lemma [7751 which follows 
this proof, i X] 't[k,m)b^j'^^ \ is a basis for the trivial component of Res^ . Ind^i?. 

We can complete it to a basis of Ind^i? by adding any bases for the non-trivial irreducible 
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components of Res^ Ind^i?. The exact form of this completion does not affect the 
quotient ^/ind%R (nor the transplantation) since the functions that correspond to these 
basis elements vanish on the edge e\ We denote this basis of Indji^i? by 03*, and this 
finishes the description of the bases and the representatives by which the two quotients 
are constructed. We summarize them and the encoding\decoding relations that they 
imply: 

• ^^R'. We have je*^*'^-* = t^^^e*}^^^^^ as representatives for the orbits '^^'/h. The 
resulting edges of the quotient are jej*'^''! ^^^^^ , where for each fc, the edges 



i<j<dt 



\ e^*'*^^ \ correspond to the basis \ 6^*'^^ \ . Let / be a function on ^/r. 



For each fc, we can decode its restrictions to the edges < e - 



into the 



l<j<dk 



set of functions < /, 



l<j<d 



whose restrictions to He^'''''^ transform according 



to the representation R given by the basis The decoding is given by 




(i,k) 



l<J<dk 



dk< 3 <d 

• ■'^/ind^-R! Obviously, we have representative for the orbit ^^'/g. 

The resulting edges of the quotient are ; and they correspond to 



l<3<df. 



the basis { ^ Yj t{k,m)b'j 

m=l 



l<k<l 
l<j<dfe 



restrictions into the edges ^^kj-^ 

1 ac 

l<J<dk 



For {p, a function on ^/ind'§R, we decode its 
,^ by the set of functions < > ^^^^^ whose 



i<fc<i 

l<j<df, " ' l<j<d 

restrictions to Ge* transform according to the basis ^\ The decoding is given by 



dk < j < d 

We recall that we wish to find an isomorphism between $r/ji(A) and $r/i„dgfl('^) for every 
A G C. In order to do that we compose the "^/ij-decoding" with the "r/indgi?-encoding", 
but in the middle we need to establish a bijection between the sets of functions whose 
restrictions to Ge* transform according to IJi=i those whose restrictions to 
Ge^ transform according to 03*. The bijection suggested by the bases we have chosen is: 



1 



t{k,m)fj 



\i,k) 



(7.13) 



m=l 



By (17121) . VI < m < n^, tV^.e' = t'^t = e^^'^\ so that fl7:T3D simplifies to 



'^{k,3) 



^{k,m) 

■n-k 



1 V^. 7{i,k) 



m=l 



Uk 

rik ^ ^ 



Xi,k) 



i{i,k) 



(7.14) 
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Having set the stage, the transplantation is simply given by: 



e'- ■' 



g(i.k) ^ 



J 

□ 



IS a 



Lemma 7.5. Following the notations of the proof above, ^(fc,m)^(fci) 

i<j<dfe 

basis for the trivial component o/ReSg Ind^i?. 
Proof. Let 1 <k <l and 1 < m < n^. 

Gsit(k,m) = [J l^gi t(fc^m) |^^(fc,m')-^} ~ LJ ^C^'™') {^(fc^m') ^C^^''") fl 

m'=l m'=l 

= U t(^k,m') {tlk,ni) Gs, t(^k,m)[^H^ = [j t(fc,m') jCi-l^^^-, Q } 
m'=l m'=l 

= [J t{k,m')H^~i^^^i = [J t(^k,m')H^(i,k). (7-16) 
m'=l m'=l 

The first equality is due to G^i ti^k,m)H = IJm'=i t{k,m')H and the third one is by fl7.12p . 
In particular we observe that for every k, \G^i\ = nk\Hg(i,k)\. For fej*''^^ with 1 < j < d^, 
we have by fl7.16p 

' ' geG-i ' ' m'=l /le/i'.f^.fc) m'=l 

where the last equality is since H^(t.k) acts trivially on b^^''^\ Obviously, the l.h.s of (17.171) 



t{k,m)b^j'^^ \ i<7n<4 is a basis for Ind^-R. To finish the proof we show that its size, 

l<j<d 

Y^k=i '^k-, equals (Res^ .Ind^-R, 1 



is invariant under the action of G^i. Therefore, the set < ^ ^ t{k,m')bj \ lies 

l<i<dfe 

in the trivial component of Res^ Tnd^i?. In addition, it is linearly independent since 

i<fc: 

<77l< 

i<j; 

^Res^ Ind^i?, l\ . 

If i? is a representation of H < G, then for any g E G we also have R 
as a representation of gHg^^, with ghg'^ G gHg^^ acting as h E H . Mackey's 
Decomposition Theorem ([38], section lOB) uses this to relate induction and restriction 
in the following way: 
I 

Resg Ind|i?^ ^Indf^; ,_iRes!'f'=' .R, 

k=l 



lilt ^ 

= 1 

which shows that 



Resg_Tnd|i?, l) =J2 ( IndJ^ t-^^^'?'" t-^' ^ 

^ k=l 



^ \ ''^!s(i,k)^k ^kJ^-(i,k)T^k / r< 
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□ 



Corollary 7.6. // G acts on F anc? Hi,H2 are subgroups of G with corresponding 
representations Ri,R2 such that Ind^^i?i = Ind^^-^a? then ^/ri and ^/r2 which are 
constructed by the recipe described are transplantable. 

We end this section by demonstrating the process of finding the transplantation 
for our basic isospectral example (figure [2]). The representations (13. ip . (13. 2p that were 
used to form those quotients are one dimensional and we denote their bases by i? = {b}, 
03 = {b}. We choose the representatives {e, cr} to form the induction Ind^*i?i whose 
basis we choose to be {eb, ab}. 

The matrix representation corresponding to this basis is 



io\ f 1 \ 2 /-io\ 3 /o-i 

1' ^^-10' ^ ^ -1 h 1 



-10\ /0-l\ 2 /10\ 3 (01 



.(7.18) 



We choose the same representatives {e, a} and get lnd^'^^R2 with the basis {eb, ab} and 
the following matrix representation 



10\ /01\ 2 /-10\ 3 (0-1 

1' ^ -1 ' ^ ^ -1 ' ^ ^1 



01\ /-10\ 2 /0-l\ 3 (10 

1 ' ^ 1 ' -1 ' -1 



(7.19) 



The relation between these two bases is 

eb = Vv^(eb - ab) 
ab = Vv^(eb + ab) . 

Let A G C. Let / G $r/fl.^(A), ip G ^r/ji^{X). Decoding these functions we get 
/i G $p^(A), ipi G $p^(A). Note that in the general construction we also have a 
superscript index denoting the number of the edge on which the function is defined. 
However, in this example we have a free action on the edges so that we can choose a 
single basis {B for ^/ri and *B for ^/r2) to work with. Using the general notation we 
therefore have VI < z < 6, fi = fl and similarly for ^i. Forming the inductions 

we have < e/i, afi ^ C $p ^ (A), {ev^i, a^pi) C $p ^ (A). These functions 
correspond to the bases {e6, ah} , {eb, ab} and transform accordingly. We now describe 
the transplantation by expressing / in terms of (p (the edges' labeling is given in figure 
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t » t •/) 



(b) 



Figure 10: (a) The graph ^/Ri with its edges labeled, (b) The graph T whose relevant 
edges are denoted, (c) The graph ^/R2. 
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The transplantation on all the other edges is obtained similarly and the result is 
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8. A gallery of isospectral graphs 

We present three interesting examples of isospectral graphs. The first example shows 
that once the algebraic conditions stated in theorem 14.31 or in corollary 14. 41 are satisfied, 
we can choose any graph that obeys the group symmetries and obtain isospectral 
quotients. The second example shows a symmetry group larger then 1)4, as well as 
the effect of using different choices of representatives while constructing the quotient. 
The third example demonstrates the special case of a group that acts freely on a graph 
which enables us to get quotients whose boundary conditions are all of the Neumann 
type. 

8.1. G = D4 — Another Ghoice for T 

We present another example using the group D4, to emphasize that it is possible to 
obtain isospectral quotient graphs from any graph that obeys the symmetry of D4. We 
now denote the Cayley graph of with respect to the generating set S = {a, r} by 
r (figure [H]). Since a Cayley graph can be constructed for any discrete group, this 
example can easily be generalized to create other isospectral sets. The set of edges of 
r is then E = {{g, gs)}g(zD4^^s£S, and the action of D4 on the edges of T is by Wh G 1^4, 
h{g, gs) = {kg, hgs), where s E S. Taking the same subgroups of D4 and representations 
as in sections [3] and O we obtain three isospectral quotient graphs ^/Ri, V^2, V-^a, shown 
in figure [121 The representatives chosen to construct the graphs are: 



v/H, = {e,a} 
v/H, = {e,cx} 
v/hs = {e, r} 



{(e,(T), ((T3,e), (e,r), {(T,Ta^)} 
{(e,a), {(T^,e), (e,r), {a,Ta^)} 
{(e,(T), (rc^^r), (r,e), (e,r)} 




Figure 11: The Cayley graph of D4 with the generating set {o", r}. The lengths of some 
edges of this quantum graph are marked. 
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(b) 



(4) (3) 
(3) (4) 

(2) (1) 




A 



A 



V2 










-i 1 1 1 



(c) 



Figure 12: The three quotient graphs obtained from the Cayley graph, (a) ^/Ri-, {hY/R2, 
(c) r/flg. The matrices describe the boundary conditions at the various vertices. The 
diagrams on the left relate the edges to the rows of the corresponding matrices, numbered 
1-4. 



8.2. G — Oh — Freedom to Ghoose Representatives 

We now consider a larger group, and treat Oh, the group of symmetries of the cube. 
The group Oh contains forty-eight elements: the identity, twenty-three rotations, and 
twenty-four reflections. A fundamental domain for the action of Oh on the cube is a 
tetrahedron. Each of the forty-eight tetrahedra that together form the cube has one of 
its faces along the external boundary of the cube, and its other three faces inside the 
cube. We begin by describing a graph, denoted by F, that obeys the symmetries of Oh- 
We consider star graphs with three edges, such that each tetrahedron contains one such 
star graph, and the edges of the star graph are connected to the centre of the interior 
faces of the tetrahedron (one edge goes to each face). This star graph is constructed 
to have three edges of different lengths, labeled a, 6, and c. Any two star graphs in 
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Figure 13: The isospectral graphs obtained from the cube, (a)-(c) are all ^/Ri, formed 
by choosing different representatives for ^/o and (d) is the graph ^/R2- Unless stated 
otherwise, Neumann boundary conditions are assumed. The boundary conditions of the 
marked vertices are given by the matrices on the right. 
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neighbouring tetrahedra are then connected at the center of their common face. T is 
the union of these star graphs. It is a three-regular, bipartite graph, with forty-eight 
vertices, and seventy-two edges of lengths 2a, 2b and 2c, each with equal multiplicity. 
We will use the subgroups 0,Td < Oh, along with appropriate representations, to form 
the isospectral quotient graphs. 

We first consider the subgroup O < O^, known as the octahedral group, which 
contains the identity and the twenty-three rotation elements. We take representatives 
for ^/o and ^/o. Since O is a subgroup of index two, one possiblitiy is to choose the edges 
and vertices contained in two tetrahedra as representatives. O has two 1-dimensional, 
one 2-dimensional, and two 3-dimensional irreducible representations [36]. We work 
with the two dimensional representation, using the basis given in [36] and we denote it 
as Ri. The quotient graph ^/ri is shown in figure fTSTa). As noted in sections El 17.41 it 
is possible to make different choices for the representatives. We present three different 
quotient graphs for ^/ri, corresponding to three different choices of representatives, in 
figure fT3ra)-(c). 

We now examine the subgroup T^ < Oh- The vertices of the cube consist of two 
sets, each of which forms the vertices of an equilateral tetrahedron. contains all the 
elements of Oh whose action does not mix between the two sets. It should be noted 
that O = Td = S4. In particular, and O have the same irreducible representations. 
We will take the matrix representation that was used for O, and use it for T^, denoting 
it by i?2- The quotient graph ^/r2 is shown in figure [T3](d). 

The isospectrality of the quotient graphs formed from O and Td follows from the fact 
that Indg'^i?! = Ind^^'^i?2- 

8.3. G = — A free action on T 

We now consider an example that demonstrates the possibility of constructing quotient 
graphs that have only Neumann boundary conditions. We show that this can be achieved 
even when dividing by multidimensional representations. We consider the graph F, 
shown in figure [141 which is symmetric under the action of D-^ = S3. The group 
acts freely on both the vertices and edges of F. This is ensured by having the circles 
perpendicular to the plane of the triangle. We have D3 = {e, a, a^, r, ra, ra^}, where a 
is a rotation of 2^/3 about the axis perpendicular to the plane of the triangle, and r is 
a rotation of vr about the height of the triangle. 

Rather than working with representations of various subgroups, we use various 
bases of a representation of the entire group to create the quotient graphs. D3 has 
three irreducible representations: the trivial representation Si, the sign representation 
5*2, and the standard representation 5*3, which is of dimension two. We take the eight- 
dimensional representation R = 3Si + S2 + 2S3. We begin by choosing a basis such 
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Figure 14: Two visualizations of the graph T which is symmetric under D3. (a) A three 
dimensional visualization. The loops are perpendicular to the plane of the triangle, a 
is a rotation of '^'^/3 about an axis perpendicular to the plane of the triangle, and r is a 
rotation of vr about a height of the triangle, (b) A two dimensional visualization whose 
edges' lengths are marked. 



that the matrix representation of R is block diagonalized as follows: 



\ 



Si 



Si j 

This choice of basis is not unique. We therefore specify it, by identifying the top block as 
the regular representation of D3, and taking the standard basis for this representation. 
Due to the fact that the matrix representation is block diagonalized, the quotient graph 
created consists of three disjoint graphs. Each graph corresponds to one block, and 
these are shown in figure [T5l 



Note that this representation consists of permutation matrices. This, together with the 
fact that the action of on F is free, ensures Neumann boundary conditions on the 
quotient graph. We now choose a new basis, such that the matrix representation of R 
takes the form: 



/ Si + S, 



\ 



Si + S-i 



\ 



Si + 5*2 



Again we must further specify the choice of basis. Each of the first two blocks is the 
permutation representation of the symmetric group of size three, and we choose its 
standard basis. For the third block, we choose the basis in which: 




r 1-^ 
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(a) 



(b) 



Figure 15: The first quotient graph ^/Si+S2+2S3 U^/Si U^/Si. All vertices have Neumann 
boundary conditions, (a) The graph formed from the first block of the representation. 
Note that in this case, we have recovered the original graph, (b)-(c) Two copies, coming 
from the two blocks containing the trivial representation. 



Using this basis, we again obtain a quotient graph consisting of three disjoint graphs, 
corresponding to the three blocks of the matrix representation, as shown in figure [I6l 
This quotient also has only Neumann boundary conditions for the same reason. The 
two quotient graphs, namely ^/si+S2+2Ss U ^/si U ^/si and ^/si+Ss U ^/si+Ss U ^/si+S2, are 
isospectral. 



(a) 




(b) 

Figure 16: The second quotient graph ^/Si+S^ U ^/Si+Ss U ^/Si+S2. All vertices have 
Neumann boundary conditions, (a)-(b) The graphs formed from the first two blocks of 
the representation, (c) The graph obtained from the third block. 

Remark. As a matter of fact, all the quotients obtained in this subsection can also be 
obtained as quotients by the trivial representations of D^s subgroups, as 



ind^n 



^1 
Si 



So 
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9. Drums and Manifolds 

We now apply theorem (14. 3p to other objects. In particular, we reconstruct some existing 
examples using our method, and comment on some new results. 

9.1. Jakobson, Levitin, et al. 

We begin by examining the isospectral domains presented by Jakobson et al. and Levitin 
et al. P [ID]- It is possible to recover all the isospectral examples described in these 
papers as quotients by representations with isomorphic inductions. We consider first an 
interesting example, consisting of the four isospectral domains shown in figure [I7| (See 
also figure 7 in p!0]). 




'^/Ri^R2 '^/R2(S>Ri V^a^-Ra 

Figure 17: The four isospectral domains presented in [TOj, obtained as quotients of the 
torus T (figure [T8]) . solid lines indicate Dirichlet boundary conditions, and dotted ones 
Neumann. 

We take a torus, T, with sixteen diamond-shaped holes, as shown in figure [181 We 
can express the torus as x T^, where T^, Tf, are circles with circumferences a and b 
respectively. Taking a rigid action of D4 on each of the circles and T^, we obtain an 
action of G = -D4 x D4 on the torus Ta x T^: the action of the first D4 in the direct 
product is by horizontal translations and reflections, and that of the second is by 
vertical ones. For example, the element [a, e) transforms T rigidly onto itself so that 
the four columns of diamonds are cyclically shifted by one, and the element (r, e) swaps 
the flrst column with the fourth, and the second with the third. Similarly, the action of 
elements of the form (e, g) is by transformations that permute the rows of diamonds. 

Adopting the notations of section [3l we examine the subgroups {Hi x Hj}.j^^-^^^y of 
D4 X D4 and their corresponding representations {Ri ® Rj}. j^^-^ Since Ind^^Ri = 
Ind^^-^25 cill of these representations have isomorphic inductions in D4 x D4. Applying 
corollary 14.41 we obtain that {V-^^'^^ilj jG{i 2} isospectral (flgure [TTjl . 
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Figure 18: The torus T, of which the four domains in figure [T71 are quotients. 
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Figure 19: (a) The information we have on / G $^i®^^(A). The funct ion vanishes along 
the sohd hues and its normal derivatives at the dotted lines are zero, (b) The quotient 
planar domain '^/Ri®R2 which encodes this information. Solid lines stand for Dirichlet 
boundary conditions and dotted lines for Neumann. 



To demonstrate how the technicalities work for these domains, we construct Yi?i®-R2. 
From and (D, we obtain information on / e ^^^'^^^{X): we see that / vanishes 
on the axes of reflection of (r, e) , as it is anti-symmetric with respect to this reflection, 
and similarly for (e,rcr). Since / is symmetric with respect to the axes of (rcr^,e) and 
(e,r(T^), its normal derivatives with respect to these axes are zero. This is summarized 
in figure [TOT a) and the resulting quotient is shown in figure [TOTb). Notice that the two 
parts of figure [19] serve the same purpose as those of figure [3l i.e., to demonstrate how 
the information on a function belonging to a certain isotypic component is encoded 
in the boundary conditions of the quotient. We end by remarking that all the 
constructions demonstrated in section [5] can be applied analogously to T to enlarge 
the isospectral quartet in figure [T71 However, the quotients obtained from most choices 
of representations and bases will not be planar domains, or even manifolds. For example, 
consider the representation of H^, given in (15.11) . The quotient '^/ri^r-j, (figure [20!) is 
a cylinder with Dirichlet and Neumann conditions at its boundaries, and a "factor of i" 
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condition along a section line normal to the boundaries (compare with the quotient ^/rs 
introduced in section [5]). This quotient is a manifold with a singularity. Other types 
of singularities may arise when considering quotients with respect to multidimensional 
representations. 



D 




Figure 20: The quotient '^/Ri(^R3. SoHd lines stand for Dirichlet boundary conditions 
and dotted lines for Neumann. There is a factor of i between the values and derivatives 
at the top and bottom edges. 



We now examine a more general example presented in [TO] . Consider the cylindrical 
drum D, shown in figure [211 In the figure, the left and the right edges are identified. 
The three dots imply that the basic pattern which appears in figure [22] is repeated and 
D consists of 8n copies of it. © is symmetric under the action of the dihedral group -D4„, 
where a rotates the cylinder and r is a reflection whose axis is shown in figure [211 



A 



rcr 



TCT 



_4h-1 



TCT 



TCT 



.4n-2 



TCT 



Figure 21: The cylindrical drum D. The left and right edges are identified. Solid lines 
stand for Dirichlet boundary conditions and dotted ones for Neumann. The reflection 
elements in D^^n and their indicated. 



I i 



Figure 22: A fundamental domain of U) under the action of D^^n. 
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(b) - »/a 



Figure 23: (a) The information we have on / E $p^(A). The function vanishes along the 
solid lines, and its normal derivatives at the dotted lines are zero. (b),(c) The quotient 
planar domains ^/Ri and ^/R2- Solid lines stand for Dirichlet boundary conditions and 
dotted ones for Neumann. These domains illustrate the case n = 2. 



We consider the subgroups Hi, H2 < D^n, where 

H, = {e, r, Ta^\ a^^} 

H2 = {e, ra, ra^-+\ a^"}, 
equipped with the one- dimensional representations 

i^l:|e^^(l), r^(-l), ra^" h-^ (1) , a^" h-^ (-1) | 

i?2 : { e ^^ (1) , TG^ (1) , ra^^+i ^ (-1) , ^ | _ 

We find that Ind^*"i?i = Ind2^"i?2, and therefore we can form two isospectral quotients, 
each one a quarter the width of the original drum, with Dirichlet boundary conditions 
on one side, and Neumann on the other. In figure [23] the drum D and the resulting 
quotients are shown for the case n = 2. In conclusion, we have provided an alternate 
proof for theorem 4.2 in [TUj. In fact, this proof is valid for any number n, whereas the 
original theorem in [10] treated only the case n = 2^ . The reader might wonder why 
did we limit our attention to D^^- The answer is that if we consider the more general 
Dn, then in order to define Hi and H2-, n must be even. If n is even, but not a multiple 
of four, then the proof works, but in that case Hi and H2 are also conjugate in D4, so 
that the isospectral domains thus obtained are also isometric. 

9.2. Gordon, Webb and Wolpert 

The famous domains of Gordon, Webb and Wolpert, originally presented in [6l [7], can 
be similarly constructed by our method. Buser et al. [8] have shown that they can 
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be constructed as quotients of the hyperbohc plane. This is done by considering an 
epimorphism from the symmetry group of this plane onto PSL3 (F2), and taking the 
inverse images of two subgroups in PSL3 (F2), each isomorphic to 5*4. In our formalism, 
the drums of Gordon et al. are obtained (with Neumann boundary conditions) from 
the quotient of the hyperbolic plane by the corresponding trivial representations of 5*4. 
Using the sign representation of S4 instead we obtain the same drums with different 
boundary conditions (figure [2^ . The conditions of corollary 14.41 are satisfied in this 
case as well, so that this pair of drums is isospectral. A more detailed explanation is 
given in [27] . 




Figure 24: The isospectral drums of Gordon et al. with new boundary conditions. 



9.3. Chapman's two piece band 

In there are no Sunada pairs, i.e., there are no nonconjugate subgroups A,B < D4 
whose trivial representations satisfy Ind^"!^ = Ind^'^l^. However, by simple arithmetic 
we can extract isospectrality even from these "basic ingredients" . 

n 

We note that if {Ri}^^^ are representations of G, then the disjoint union |J t/r. is 

1=1 

n 

isospectral to ^/ ^Rt, as we have $ n (A) = ^r/n^ (A), so that 

i = l U ^/^i i = l 

r/n^ = H {XR.,X'l>rW)G = (^eR^'^*r(A) ) = ^^r/ W 
1=1 ' i=l \ 1=1 ' / G 

(this was implicitly manifested in section [8731) . Combining this with theorem 14.31 we 
obtain the following: if {Hi},\^Hj^ are finite sets of subgroups of G, with corresponding 
representations {Ri}, {R'j}^ such that 0^1nd^.i?j = 0^- Ind^/i?^-, then [J^^/Rt is 
isospectral to [Jj^/R'j- 

We recall the subgroups Hi,H2,H^ < D4 from sections [3], [5] and consider also 
H4 = {e, ro"^} and = {e, rcr}. Even though no two inductions among {ind^'^ln,} 
are isomorphic, we do have that Ind^^^lni © I^d^^lf^g = Indj^^l/^j © Ind^^^l/^^. By 
the observation above, if D4 acts on some object F, then V^^i U V^^s is isospectral to 
r/i,,ur/i,^. 
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Figure 25: (a) A square S, on which acts, with the axes of reflection elements marked, 
(b) Fundamental domains for the actions of Hi = (^T,Ta'^'j, H2 = (ro", ro"^^, = (^ro"^^, 
H5 = {tu) on §. 



In figure [231 part (a) displays a square S, of side length 2, on which acts, and 
part (b) displays fundamental domains for the actions of D^s subgroups {Hi}^^^ 245- If 
we supply § with Neumann boundary conditions, we obtain that {^/ih,} are the domains 
shown in figure [251(b). also with Neumann boundary conditions. The isospectrality of 
^/ihj U ^/ih5 and ^/ihj U ^/ih4 is a known example constructed by Chapman |29j, which 
shows that "one cannot hear the shape of a two-piece band" Chapman obtains 

this isospectral example by manipulating the drums of Gordon et al. — he enlarges the 
number of connected components by cutting the basic building block, and then cancels 
out identical components in the two shapes. Chapman shows that these domains are 
also isospectral if the Neumann conditions at the boundaries are replaced by Dirichlet 
ones, and we would like to establish this as well. A reasonable guess would be to try D, 
a square identical to S but with Dirichlet boundary conditions. This, however, leads to 
isospectrality of domains with mixed boundary conditions (figure [2^ . 



yi 



Hi 



Figure 26: Various quotients of O, the square in figure [251(a) with Dirichlet boundary 
(Solid lines stand for Dirichlet boundary conditions and dashed ones for Neumann). 
U ici/iHg is isospectral to 1^/1^2 U 10/1^4. 



In figure [261 the Dirichlet edges are the remnants of D's boundary, whereas the 
Neumann edges are the traces of the reflections by which the quotients were taken. One 
is thus led to consider the following representations: 

^'i:{e^(l), rH^(-l), ra^ (-1) , ^ (1) } 
K--{e^{l). Ta^{-1), ra3^(-l), ^ {I) ] 
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K: {e^{l), ra' ^ (-1) } 

K: {e^ (1) , ra ^ (-1) } . 

In {R'j}, the reflection elements are sent to (—1), so that the corresponding quotients 
encode functions which are antisymmetric at the corresponding axes, hence vanishing 
along them. The quotients {"^/fi-jj^i 2 4 5 indeed the domains presented in figure 
1251(b) with Dirichlet boundary conditions, and upon verifying that lnd^lR[®lnd^lR'^ = 
Ind^^i?2 © Ind2;^i?4, we obtain the isospectrality we have sought. Note that by taking 
{^/^i}i=i 2 4 5' would have again obtained the mixed isospectral example in figure [26| 
but with the Dirichlet and Neumann conditions swapped. We end by remarking that as 
usual, one can extend the isospectral families above by considering other bases for the 
representations Ind^^l//j © Ind^*!//^ and lnd^'^R[ © Ind^^i^g, but that in the general 
case the objects thus obtained will no longer be planar domains. 

10. Summary and open questions 

This paper describes a method which enables one to construct isospectral objects, such 
as quantum graphs and drums. The algebraic component of the underlying theory 
suggests theorem 14.31 and corollary 14.41 as our main tools for producing isospectral 
objects. Their assumptions are less strict than those of Sunada's theorem [2] and this 
allows more degrees of freedom in the isospectral search. Another ingredient is the 
assembly process of the so called quotient graphs, whose construction accounts for yet 
more liberty. We found that for a graph F with a symmetry group G and a representation 
R of the group, we have a variety of choices to make for the fundamental domain of the 
action of G on F and also for the basis with respect to which R is presented. These 
different possibihties yield (possibly) different quotient graphs ^/r, all isospecral to each 
other. We wish to offer two perspectives on this dizzying freedom. On the one hand, it 
invites us to test the strength of the method. Namely, given two isospectral objects, can 
it be decided whether they arise as quotients of some common object? We touched this 
question so far merely by reconstructing some known isospectral objects in terms of our 
method. On the other hand, it prompts one to classify these sources of isospectrality 
and understand the interrelationships between them. A fundamental question in this 
context is whether for R as above, there exists a choice of basis for Ind^R which makes 
^/ind^R the same as ^/r (rather than just isospectral). We saw a demonstration of this 
in section [5] where ^/ri was the same as ^/indg^iJi with an appropriate choice of basis. 
If this is always the case, it means that the ability to change between different bases 
is a more fundamental source of isospectrality. This should not cause one to abandon 
theorem 14.31 and corollary 14. 4[ Their role in such a case would be to indicate favorable 
bases for the construction. First, they make the practical assembly of the quotient easier 
by offering lower-dimensional representations to divide by. Secondly, the quotient of a 
manifold by a multidimensional representation is seldom a manifold, so one is led to seek 
one-dimensional representations with isomorphic inductions (see examples in section [9]). 
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Pondering over the quotient graph, which stands at the heart of our method, 
we are led to inquire how its various properties are determined by the construction. 
Among these are topological ones, such as the number of connected components and the 
number of independent cycles of the graph. Others relate to the nature of the boundary 
conditions, which in turn determine the qualities of the differential operator on the 
graph. Of specific importance are conditions which guarantee that the quotient graph 
has only Neumann boundary conditions, or alternatively, boundary conditions which 
ensure the self-adjointness of its Laplacian. This issue was addressed in proposition 17.11 
and the example in section 18.31 but still awaits further investigation. 

Quantum graphs are the focus of this paper and they obtain a thorough treatment. 
One reason for this, which was already mentioned, is that under fairly mild conditions, 
the resulting quotient object is also a quantum graph. The other reason is that it is 
relatively simple to give a rigorous description of their construction (see section [6]). 
However, we have demonstrated in section [9] that the method is also applicable to 
manifolds and drums, and it is desirable to examine the possibility of obtaining other 
isospectral objects as well, e.g., discrete graphs. 

Another interesting application would be to relate the construction method to the 
spectral trace formula for quantum graphs. Specifically, we would like to show the 
equality of the trace formulae for isospectral graphs just by examining the way in which 
they are constructed. This includes a comparison of the total lengths of the graphs and 
the lengths of their periodic orbits. A similar question for isospectral planar domains 
is discussed in [30]. We propose to treat this issue by returning to the origin of the 
spectral trace formula for quantum graphs [HI [18], which was developed by describing 
the boundary conditions in terms of scattering matrices. Therefore, it may be worthwhile 
to work out an isospectral theory, analogous to the one described in this paper, but 
stated in terms of scattering matrices. Such an approach may also pave the way for a 
similar isospectral theory for discrete graphs, as a spectral trace formula for them was 
recently developed using scattering matrices [3T] . 

We end by returning to Kac's question and asking what can one do when hearing 
the shape of a graph (drum) is not possible. One answer concerns the field of counting 
nodal domains of the Laplacian's eigenfunctions. Some new works investigate the ability 
to resolve the isospectrality of discrete graphs, quantum graphs and various manifolds 
by counting their nodal domains [231 1321 [331 ^M- specific method of doing so by 
relating the nodal count of an isospectral pair to its transplantation was developed in 
[26] . The theory presented in this paper and the transplantation it yields can perhaps 
lead to a general method of isospectrality resolution. It may therefore be further asked 
whether one can count whatever cannot be heard. 
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Appendix A. A short review of required elements of representation theory 

Let G be a finite group. A d-dimensional representation of G, denoted by R, consists of 
a vector space of dimension d equipped with an action of G, which is described by 
a homomorphism pR : G ^ GL{V^), i.e., Vg'i, g2 G G, PR{gig2) = PR{gi)pR{g2)- Once 
a basis for is chosen, one can also think of pr as a homomorphism into GL^ (C). 
Pr is called the structure homomorphism of i?, and the vector space the carrier 
space of the representation. We alternatingly use R, pR and when refering to the 
representation. 

The character of a representation R is defined as : G — C, Xnid) = ^^(P-R(fl'))- 
We will also use the notation xv for the character of a representation R whose carrier 
space is V. There is an inner product defined on the characters by 

(Xiii, Xii2>G — -^^^XrA9)XR2{9)- 

A representation R is called reducible if there exists a nontrivial subspace of the carrier 
space which is invariant under the action of the group. Otherwise it is irreducible. Up 
to isomorphism, any finite group G has a finite number r of irreducible representations, 
{'5'i}[=i- We often use S to denote irreducible representations, and R for general 
ones. The characters of the irreducible representation obey orthogonality relations 
Vi, J e {l..r} , (xsi, XSj)g = Sij. 

Two important notions that are used throughout the paper are the restriction and 
the induction of a representation. Let if be a subgroup of a group G. Let R be 
a representation of G. Then the restriction of R from G to H, denoted Res^-R, is 



. In particular, the dimension of Res^i? 



described by \/Rcs^r _ yR PResffi — PR 
is equal to that of R. 

Now, let i? be a representation of H. We describe the induction of R from H to G, 
denoted Ind^i? . We start with a carrier space W — for the representation R of 
H and construct a vector space V — V^'^^"^ which carries a representation of G. We 
choose representatives for the left cosets of in G: {ti}^^-^, where n — [G : H]. For each 
i e we form a space denoted tiW, which is isomorphic to W. The prefix ti of the 
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elements tiV G tiW is currently only an abstract notation without an actual meaning. 
However, as one may expect, the vector tiV will obtain, via the described construction, 
the meaning of the action of tj G G on f G V^. The desired carrier space V is defined 
to be = 0"=]^ tiW and we now equip it with an appropriate action of G. Let v E V, 
g E G. There exist unique {wj}"^^ C W such that v = '^i'^iy ^^"^ there exist unique 
hi E H and a (z) G {l..n} (which depend on g) such that gti = t„(^i)hi for i G {i..n}. We 
are motivated by the "identity" gv = Yl^=i Q'^i'^i = SILi '^a-{i)hiVi to define 

n 

g ■ V = '^ta{i)ihiVi) . (A.l) 
i=l 

Where hiVi is computed according to the representation R from which we started. Note 
that the action of g permutes the subspaces {tiW}^^^ among themselves and in addition 
manifests the action of a corresponding H element (determined by g and ti) on each 
such subspace. It is left for the reader to check that this indeed defines an action, i.e., 
^gi,g2 G G, ((71(72)"^ = 5'i(5'2'y), and that the obtained representation does not depend 
(up to isomorphism) on the choice of representatives. Also, V contains a subspace 
which as a representation of H is isomorphic to R, namely, the tiW which corresponds 
to the trivial if-coset. In contrast to the preservation of dimension of the restricted 
representation, we have dimlnd^i? = dimi? ■ [G : H]. 

It is known that the character completely identifies a representation and it 
is therefore useful to describe the characters of the restricted and the induced 



, and explicit calculation gives 

H 



representations. Obviously, Xrcs^b. = Xr 

n 

Xindg/?(^) = XnitJ^gU). (A.2) 

i=l 

where XRi^7^9U) = for all t^^gti ^ H. 

Let us demonstrate this by referring to section [3] and calculating the induction of 
the representation Ri from Hi to G = D4. We choose the representatives {e, a} for the 
cosets of Hi in D4. The character of the induced representation is therefore 

Xindg^ iJi (e) = XRiie) + XR^{e) = 2 

Xindg^iJi (f^) = XRi (c^) + XRi ((xVa) = 

Xind%^RM^) = XRii(^^) + XRM'^cr^(^) = -2 

Xindg^ rM) = XRiir) + XRr{(r^r(r) = 

Xind%^RM^) = XRiira) + XR^f^^raa) = 
In the above we have calculated the character only for chosen representatives of the 
conjugacy classes of G (which is obviously enough). We indeed obtain the character of 
the two-dimensional irreducible representation of D4 (see for example (15.71) ). 

We end by stating a fundamental theorem, known as Frobenius reciprocity (|38j. 
section lOA): Let G be a group with a representation _Ri. Let if be a subgroup of G 
and i?2 a representation of H. Then 

(XResgi?i5XiJ2)H = (XHi,XlndgR2)G- 
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This is the main element that stands behind the proof of theorem 14.31 
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